arXiv: 1505.0329Ivl [cond-mat.mtrl-sci] 13 May 2015 


Water on hexagonal boron nitride from diffusion Monte Carlo 

Yasmine S. Al-Hamdani/'^ Ming Ma,^’^ Dario Alfe/'^ 0. Anatole von Lilienfeld,"^’^ and 
Angelos Michaelides^’^’0 

Thomas Young Centre and London Centre for Nanotechnology, 

17-19 Gordon Street, London, WCIH OAH, U.K. 

Department of Chemistry, University College London, 20 Gordon Street, London, 
WCIH OAJ, U.K. 

Department of Earth Sciences, University College London, Gower Street, 

London WCIE 6BT, U.K. 

Institute of Physical Chemistry and National Center for Computational 
Design and Discovery of Novel Materials, Department of Chemistry, 

University of Basel, Klingelbergstrasse 80 CH-fOdd Basel, 

Switzerland 

Argonne Leadership Computing Eacility, Argonne National Laboratories, 

9700 S. Cass Avenue Argonne, Lemont, Illinois 60439, 

USA 

(Dated: 14 May 2015) 

Despite a recent flurry of experimental and simulation studies, an accurate estimate 
of the interaction strength of water molecules with hexagonal boron nitride is lacking. 
Here we report quantum Monte Carlo results for the adsorption of a water monomer 
on a periodic hexagonal boron nitride sheet, which yield a water monomer interaction 
energy of —84 ± 5 meV. We use the results to evaluate the performance of several 
widely used density functional theory (DFT) exchange correlation functionals, and 
hnd that they all deviate substantially. Differences in interaction energies between 
different adsorption sites are however better reproduced by DFT. 
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Hexagonal boron nitride (h-BN) has become popular for anyone with an interest in 2- 
dimensional materials, due to a number of notable properties such as high thermal conduc¬ 
tivity, mechanical robustness and exceptional resistance to oxidation,- and not least because 
it is isostructural with graphene. Our interest has been piqued by experimental reports of 
fascinating behavior of water at h-BN such as superhydrophobicity,- water cleaning ability^ 
or generation of electric current.- These experiments have already prompted a number of 
simulation studies of water on h-BN sheets and nanostructures using both density functional 
theory (DFT) and classical molecular dynamics.-^— They have been incredibly informative 
and have helped to e.g. understand the electrical currents generated in BN nanotubes.- 
However, there is one major unknown at the very heart of any water/BN simulation 
study: we simply do not know what the interaction strength of a water molecule with 
h-BN is. DFT calculations yield a range of values for the water monomer adsorption en¬ 
ergy depending on the exchange-correlation (xc) functional used^^- and force fields rely on 
interaction parameters htted to particular xc functionals or to experimental data such as 
contact angles for macroscopic water droplets.-"— If htting to experiment one needs to be 
certain that the experimental conditions are exactly known; recent lessons learned for wa¬ 
ter droplets on graphene reveal that contact angle measurements are incredibly sensitive to 


surface preparation conditions and levels of cleanliness.— — 

The lack of well-dehned reference data for water on h-BN is representative of a much 
broader problem: there are very few systems for which accurate water monomer adsorp¬ 
tion energies have been established. Mainly this is because even at low temperatures water 
molecules cluster into larger aggregates making the determination of monomer adsorption 
energies with established surface science techniques such as temperature programmed des¬ 
orption or single crystal adsorption calorimetry highly challenging.—"— In the absence of 
experimental data simulations play an important role, either via explicitly correlate d q uan- 
tum chemistry approaches or quantum Monte Carlo (QMC) (see e.g. references Il9l425l). 


Indeed given recent increases in computational capacity and the fact that it can be applied 
to periodic systems, QMC has emerged as a powerful technique for obtaining interaction 
energies of molecules with surfaces^^*^ or biomolecules.—"— 

Here we report results for interaction energy curves for water on a periodic h-BN sheet 
using fixed node diffusion Monte Carlo (DMC). We obtain from this an estimate of the 
water/h-BN interaction strength of about —84 ± 5 meV at an equilibrium water-surface 
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FIG. 1. Structures of the two adsorption modes of water on h-BN considered in this study. (a,c) 
Top and side view of water above an N site of h-BN. (b,d) Top and side view of water above a B 
site of h-BN. Boron is pink, nitrogen is blue, oxygen is red and hydrogen is white. All calculations 
have been performed on periodic unit cells, with the periodic unit cell in the x,y plane indicated 
by the blue frames in (a) and (b). 


distance of ca. 3.4 A. We have also computed interaction energy curves with a variety of 
DFT xc functionals and we find that these differ significantly from DMC. Except for LDA, 
of the functionals considered those that do not account for van der Waals underbind and 
those that do, overbind. DFT based predictions of the equilibrium adsorption height are 
much better with several functionals coming within 0.2 A of DMC. In addition, based on 
DMC and DFT calculations of water on h-BN in a second metastable adsorption structure, 
we find that many of the xc functionals do reasonably well in predicting the relative energy 
difference between the stable and metastable adsorption structures. 

Two different levels of theory have been used in this study, fixed node DMC and DFT. A 
standard computational setup has been used for each and so we only discuss the key features 
here. 


QMC calculations were undertaken using the CASINO code,— with Slater-Jastrow type 
trial wavefunctions in which the Jastrow factor contains electron-nucleus, electron-electron, 
and electron-electron-nucleus terms. We used Trail and Needs pseudopotentials^i^ for all 
atoms, in which the Is electrons of B, N, and O were treated as core. This set-up for the DMC 


calculations is similar to the one used in reference 


28 


where water adsorption was examined 


on 1,2-azaborine and agreement between DMC and coupled cluster with single, double and 
perturbative triple excitations (CCSD(T)) to within 9 meV was obtained. The initial single 
particle wavefunctions for use in DMC were obtained from DFT plane-wave calculations 
using the PWSCF package.— A standard 300 Ry energy cut-off was applied and for efficiency 
the resulting wavefunctions were expanded in terms of B-splines^ using a grid multiplicity 
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of 2.0. Trial wavefunctions were generated using the local density approximation (LDA)^ 
which has been validated for weak interactions in previous work.—>2^ After optimization of 
the trial wavefunctions in variational Monte Carlo, we used 6,553,840 walkers across 16,384 
cores for each point along the DMC interaction energy curves. The locality approximation 
was utilized^ with a time step of 0.015 a.u. which we tested against a time step of 0.005 
a.u. 

VAST 5.3.5^--^ was used for the DFT calculations, making use of projector augmented 
wave (PAW) potentials'*^ to model the core regions of atoms (again the Is electrons of 
B, N, and O were treated as core). Following careful tests, we chose a 500 eV plane- 
wave cut-off and a (4 x 4) unit cell of h-BN with 16 A between sheets, along with F-point 
sampling of reciprocal space.— The proliferation of DFT xc functionals over the last decade^ 
means that we cannot possibly consider all xc functionals or even all modern xc functionals 
designed to capture weak interactions. Rather we consider a small selection that have been 
widely used in adsorption studies. This includes the LDA, the PBE^ generalised gradient 
approximation (GGA), two hybrid functionals (PBE0^>^ and B3LYP— ^— ), and several van 
der Waals (vdW) inclusive functionals (PBE-I-D2,— PBE-I-D3,—DFT-|-vdW,— optB86b- 
vdW— and vdW-DF2^‘^). The DFT-|-vdW correction (from Tkatchenko and Scheffler) 
was applied to three xc functionals (PBE, PBEO and B3LYP). 

Results for interaction energy curves of a water monomer on h-BN in two different ad¬ 
sorption modes (Figure [T]), obtained from DMG and a range of DFT xc functionals, are 
shown in Figure [2l The interaction energy between the adsorbate and substrate is plotted 
as a function of the perpendicular distance between the oxygen atom of the water molecule 
and the h-BN sheet. The absolute interaction energy between water and the substrate were 
calculated as follows, 

= E‘f - E‘i (1) 

where is the total energy of water and h-BN at a given oxygen-surface separation 
distance, d, and is the total energy of water and h-BN at 8 A oxygen-surface distance. 
This dehnition allows the same Jastrow factor to be used for all conhgurations, including 
the reference structure in DMG. Adsorption structures were obtained from optB86b-vdW 
optimizations of water on a hxed h-BN sheet, see Supporting Information for coordinates 
of the adsorption and reference conhgurations.— We chose to use the same structures for 
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FIG. 2. (a) Interaction energy curves for water situated above the N site in h-BN as shown in the 
inset, (b) Interaction energy curves for water situated above the B site in h-BN as shown in the 
inset. The lines connecting the data points are merely there to guide the eye. 


DMC and all xc functionals because this makes for a cleaner comparison. For many of the 
xc functionals we have computed interaction energy curves with fully relaxed structures and 
the differences between the relaxed and the hxed structures are < 5 meV, except in the 
repulsive wall at short oxygen-surface separations. The hrst adsorption structure considered 
has the oxygen of the water molecule above an N site with one of the OH bonds directed 
at that N atom (Figure [I](a,c)). This is the most stable adsorption structure according to 
previous DFT studies.- The second structure has the oxygen atom of the water molecule 
above a B site with the plane of the molecule tilted away from the substrate by 128° (Figure 
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[T]^b,d)). According to our DFT calculations this is the most stable structure for water at 
the B site but ~ 20 meV less stable than the N site adsorption structure. We consider water 
adsorption at the B site to establish if the DFT site preference for this system is correct.— 

Let us now focus on the DMC interaction energy curves for water on h-BN. Because of 
the enormous computational cost of DMC we can only compute a small number of points for 
each energy curve, which limits the resolution of the curves. Nonetheless they are sufficiently 
well defined to yield an adsorption energy of —84 ± 5 meV at a height of r\j 3.4 A at the N 
site and an adsorption energy of —63 ± 5 at a height of ~ 3.2 A at the B site. The O atom 
sits slightly further away from the substrate at the N site because of the orientation of the 
molecule at this site, wherein the H atoms points to the N, forming a weak hydrogen bond 
like interaction. The relative energies of the two sites conhrms the DFT site preference but 
more importantly provides an estimate of the water monomer interaction energy that is free 
of any arbitrary choices of DFT xc functional. 

Obtaining an accurate estimate of the interaction strength between water and h-BN 
is important in its own right, however, it also provides a valuable benchmark which we 
now exploit. Here, we use our DMC derived interaction energy curves to evaluate how 
various DFT xc functionals perform for this system. Interaction energy curves from several 
functionals are included in Figure |2] and in some respects these reveal a familiar story. 
Looking at the most stable site first, LDA overbinds by predicting an adsorption energy of 
— 183 meV with the molecule 0.4 A closer to the surface than DMC. In contrast the GGA 
and the hybrid functionals underbind: PBE is ~ —45 meV, PBEO is ~ —40 meV and B3LYP 
is ~ —15 meV. The PBE and PBEO adsorption heights are fairly reliable at 3.40 A whereas 
the shallow B3LYP minimum is located at 3.55 A. More interesting are the results from the 
vdW inclusive functionals since these are in principle designed to accurately describe weak 
interaction systems. Surprisingly, we hnd that all vdW inclusive functionals considered 
signihcantly overbind this adsorption system. Specifically the adsorption energies are in 
the —140 to —170 meV range, with vdW-DF2 predicting the smallest adsorption energy 
and optB86b-vdW the largest. This overbinding also persists at large adsorbate-substrate 
distances; compare for example the DMC and vdW-inclusive DFT results at 4 — 5 A from 
the surface. The predicted height above the surface is in reasonably good agreement with 
DMC, only around 0.1 A closer to the surface for all vdW-inclusive functionals. 

Moving to the B site adsorption structure we find that systematically, with the exception 
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of PBE+D2, the interaction strength is reduced by ~ 20 — 30 nieV. This is in very good 
agreement with the DMC energy difference between these two sites (PBE-D2 predicts that 
the B site is ~ 60 meV less stable than the N site). Thus although none of the xc function¬ 
als considered come within 40% of the DMC interaction energy, the change in interaction 
energies between adsorption sites are in most cases described fairly accurately. We note that 
we have considered just two adsorption structures and considerably more work is needed to 
fully substantiate this conclusion. Moreover, at this stage we do not fully understand the 
poor performance of the vdW functionals and defer a more detailed analysis to a future 
publication in which results from on-going water adsorption calculations on BN clusters will 
also be presented. 

In summary, we have obtained DMC interaction energy curves for water on a periodic 
hexagonal sheet of BN and used these to evaluate the performance of a number of xc func¬ 
tionals. The interaction energy obtained is —84 ± 5 meV. This is clearly a small number; 
corresponding to the physisorption regime.— It is, however, about 15 meV larger than the 
value predicted by DMC for water on graphene. Interestingly many of the van der Waals 
inclusive functionals also predict a similar 15-20 meV increase upon going from graphene 
to h-BN.-i^S This suggests that although interaction energies are overestimated with these 
functionals, the relative interaction energies between the two materials are fairly well de¬ 
scribed. 
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